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ABSTRACT 

The energetics and emission mechanism of GRBs are not well understood. 
Here we demonstrate that the instantaneous peak flux or equivalent isotropic 
peak luminosity, Lisa ergs s~^, rather than the integrated fluence or equivalent 
isotropic energy, Eiso ergs, underpins the known high-energy correlations. Using 
new spectral/temporal parameters calculated for 101 bursts with redshifts from 
BATSE, BeppoSAX, HETE-II and Swift we describe a parameter space which 
characterises the apparently diverse properties of the prompt emission. We show 
that a source frame characteristic- photon-energy/peak luminosity ratio, K^, can 
be constructed which is constant within a factor of 2 for all bursts whatever their 
duration, spectrum, luminosity and the instrumentation used to detect them. 
The new parameterization embodies the Amati relation but indicates that some 
correlation between Epeak and Eiso follows as a direct mathematical inference 
from the Band function and that a simple transformation of Eiso to Liso yields a 
universal high energy correlation for GRBs. The existence of indicates that 
the mechanism responsible for the prompt emission from all GRBs is probably 
predominantly thermal. 

Subject headings: Gamma Rays: bursts — radiation mechanisms: non-thermal 
— ISM: jets and outflows 

1. Introduction 

The energetics of the central engine which powers the explosion responsible for a GRB 
are both intriguing and fundamental to our understanding of these cosmic events. The 
isotropic energy outflow at source, estimated using the integrated gamma-ray fluence, is 
enormous, up to Eiso ~ 10^^ ergs, and even if the outflow is coUimated in jets the total 
energy involved is still huge, Ej ~ 10^^ ergs. The possibility that the explosion taps a 
standard energy resevoir has been pursued by many authors following the initial suggestion 
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from Frail et al. (2001). If this total energy available were, indeed, roughly constant (or 
predictable through other means) and we could reliably estimate the coUimation, then GRBs 
could be used as a cosmological probe to very high redshifts. Bloom et al. (2003), Ghirlanda 
et al. (2004). 

Early on it was noted that, based on analysis of BATSE data, there was a correlation 
between i?p, the peak of E.F{E) where F{E) ergs cm^^ keV^^ is the observed spectrum, 
and the fiuence (Mallozzi et al. 1995, Lloyd et al. 2000). When redshifts became available 
for long bursts the isotropic energy, Eiso, could be estimated from the flucnce and the peak 
energy could be transformed into the source frame, Epz, the so-called Amati relation, a 
correlation between Eiso and Epz in the sense that more energetic bursts have a higher 
Epz, was discovered using data from BeppoSAX, (Amati et al. 2002). This correlation has 
subsequently been confirmed and extended although there remain many significant outliers, 
including all short bursts. The physical origin of the correlation may be associated with 
the emission mechanisms operating in the fireball but the theoretical details are far from 
settled (sec the discussion by Amati (2006) and references therein). More recently a tighter 
correlation between Eiso, Epz and the jet break time, tbreak, measured in the optical afterglow 
has been reported (Ghirlanda et al. 2004). This is explained in terms of a modification 
to the Amati relation in which Eiso is corrected to a true coUimated energy, E^, using an 
estimate of the coUimation angle derived from Ureak- The details of the coUimation correction 
depend on the density and density profile of the circumburst medium, Nava et al. (2006) 
and references therein. Multivariable regression analysis was performed by Liang & Zhang 
(2005) to derive a model-independent relationship, Eiso oc E^f^t'^l'^^, indicating that the 
rest-frame break time of the optical afterglow, tzbreak was indeed correlated with the prompt 
emission parameters. 

Other studies have concentrated on the properties of the isotropic peak (maximum) 
luminosity, Liso ergs s~^, measured over some short time scale ~ 1 s, rather than the time 
integrated isotropic energy, Eiso- Yonetoku et al. (2004) noted a correlation between Liso 
and Epz for 16 GRBs with firm redshifts. A correlation between Liso and the spectral lag 
was first identified by Norris et al. (2000) and explained in terms of the evolution of Ep^ak 
with time. The shocked material responsible for the gamma-ray emission is expected to cool 
at a rate proportional to the gamma-ray luminosity and it has been suggested that Ep^ak 
traces the cooling (Schaefer 2004). A similar correlation between Lisa and the variabiUty of 
the GRB (V) was described by Reichart et al. (2001). The origin of the Liso — V relation 
is Ukely to be related to the physics of the relativistic shocks and the bulk Lorentz factor of 
the outflow. It could be that high T outflow results in high L^so and V while lower luminosity 
and variability are expected if Toutfiow is low (see, for example, Meszaros et al. 2002). A 
rather bizzare correlation involving Liso, Epz and variability was found by Firmani et al. 
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(2006). They employed the "high signal" time, 745, as formulated by Reichart et al. (2001) 
in their study of variability, and showed that L^so oc E^f^T^^^''^^ for 19 GRBs with a spread 
much narrower than that of the Amati relation. There is currently no explanation for such 
a correlation although it may be connected with the spectral lag and variability correlations 
and the Amati relation. 

The correlation between Ei^o and Ep^ supplemented by additional empirical information 
can be used in pseudo rcdshift indicators, for example Atteia (2003), Pelangeon & Atteia 
(2006), but the intrinsic spread in the correlation and uncertainty about the underlying 
physical interpretation introduce errors, typically of a factor ~ 2. It may be possible to reduce 
the errors by simultaneous application of several independent luminosity /energy correlations, 
and extension of the Hubble Diagram to high redshifts using GRBs has been attempted, see 
for example Schaefer (2007). However, it is not clear that the correlations briefly described 
above are truly independent and there may be some underlying principle or mechanism 
which connects them all together. Recently, and more controversially, Butler ct al. (2007) 
have raised serious doubts about the validity of these correlations suggesting that it is likely 
that they are introduced by observational/instrumental bias and have nothing to do with 
the physical properties of the GRBs and hence they conclude that GRBs are probably 
useless as cosmological probes. Here we take a new look at the source frame spectral and 
temporal properties of a large number of GRBs for which we have redshifts in order to try 
and understand what really correlates with what and whether or not this can provide useful 
intrinsic information about the GRBs and what drives them. In this analysis we include 
the short-duration GRBs which may share a similar emission mechanism with long bursts 
despite probably having different progenitors. 



2. Source frame spectra of the prompt emission 

The profile of the prompt energy spectrum of all GRBs is well represented by a Band 
function (Band et al. 1993), 

B{E) = E-^^ exp{-E/E,), E < E,{I3^ - Px) 

B{E) = E-^'' exp(/3^ - fixW, - /^x)^/^"^", E > E,{p^ - (5x) (1) 

where I3x and {3^ are the spectral power law indices at low (X-ray) and high (7-ray) energies 
respectively and Ec keV is the high cut-off energy. Note that in the original formulation of 
Band et al. (1993) photon indices were used and the profile described the photon number 
density (because these are the parameters which most closely describe the detected count 
spectrum which is fitted) . Here we choose to use an energy density profile and energy spectral 
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indices. The observed total fluence is 




Ft 



tot — 



/ F{E)dE = Ntot / B{E)dE 



(2) 




ergs cm~^, where Nfot is the normahsation in ergs cm~^ keV~^ at 1 keV and Ei to E2 is the 
observed energy band. Spectral fitting of the observed count spectrum will yield values for 
j3x, P'y, Ec and Ntot- The cut-off energy, Ec, is often converted to the peak energy of the 
E.F{E) spectrum which is given by Ep = {l—f3x)Ec and the normalisation may be expressed 
as the fluence, Ftot, rather than the energy density at 1 keV, Ntot- However, the separation of 
the fluence into a normalisation term and a spectral integral is central to the development of 
the argument which follows. Table [Ugives the spectral parameters for 101 GRBs for which we 
have redshift values and a prompt light curve. The spectral parameters for bursts detected 
by BATSE, BeppoSAX, HETE-2 and Konus/WIND were taken from the references cited. 
The values for Swift bursts were derived from the BAT spectra supplemented by detections 
by INTEGRAL and Konus/WIND where available. Many of the Swift spectra (~ 40) are 
adequately fitted by a simple power law or a cut-off power law with E^. fixed. For these bursts 
a cut-off power law model was used with E^ = 150 keV (corresponding to the upper limit of 
the BAT energy band). Providing the fitted (3x < 1 the fitted function has a peak in E.F{E) 
and a value for the peak energy can then be estimated. The spectra of 7 very soft Swift 
bursts with redshifts (GRB050406, GRB050416A, GRB050824, GRB051016B, GRB060512, 
GRB060926 and GRB070419A) gave /3x > 1 and these were discarded because, for such 
spectra, we have no meaningful estimate of Ep. Such GRBs are normally designated as X- 
ray flashes (XRFs) and the exceptionally high (3x values may arise because we are actually 
observing the high energy tail {(3^) and not the lower energy power law in the Band function. 
Alternatively it may be that such soft spectra are the result of a second soft X-ray component 
which dominates in these objects. 

The equivalent isotropic energy from the source is given by 



ergs, where dL is the luminosity distance corresponding to the redshift z under some cosmol- 
ogy, Iboi{Epz, Px, P-y) is the bolometric integral of the spectral energy profile in the source 
frame, B^, taken over the wide energy band 1 keV to 10 MeV 



^T^dlNtot 



hol{Epz, Px, P-y) 



(3) 
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(1 + z)2-/3x 




(4) 



and Epz = Ep{l + z) is the peak energy in the source frame. The first term in Equation [3] is 
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the equivalent isotropic energy density, Qz ergs keV ^ at 1 keV in the source frame. 

A factor 1/(1 + zY~^^ arises because we have shifted the normahsation from 1 keV in 
observer frame to 1 keV in the source frame. The remaining factor of + z) accounts for 
the time-dilation of the duration over which the bursts are seen. It is pertinent to transform 
this to the isotropic energy density at the peak energy, Ep^ keV, in the source frame, 

Qp. = Q. exp[(i - ^x){e;^ - (6) 

ergs keV^^ so that the spectrum normalisation is specified at a characteristic energy in or 
close to the observed 7-ray energy band. We can then write Equation [3] as 

where 

E^z = exp[(/5x - l){Ep^ - l)]E^^hoi{Epz, Px, P,) (8) 

keV is a characteristic photon energy which depends on the profile of the energy spectrum 
and the limits adopted for the integration and it serves to convert from an energy density 
{Qpz ergs keV~^) at the peak of the E.Fz{E) spectrum to the total isotropic energy {Ei^o 
ergs). The isotropic energy spectrum in the source frame is given by 

FziE) = QzBz{E) = QpzeMiPx - 1){E^^ - l)]El-Bz{E) (9) 

ergs keV~^. The source frame spectra of the GRBs listed in Tabled] are shown in Figure [T] 
with the spectral energy density Qpz marked at energy Epz keV. In the majority of spectra 
the high energy spectral index is not measured but set to (3^ = !.?> which is the approximate 
average found by BATSE. Figure [2] shows the corresponding E.Fz{E) spectra in ergs. We 
assumed a cosmology with Hq = 71 km s~^ Mpc~^, A = 0.27 and Q = 0.73 to calculate the 
luminosity distance di. 



3. The Amati relation 

The Amati relation is a correlation between Epz and Eiso, first reported by Amati et 
al. (2002), and subsequently shown to be obeyed by the majority of long GRBs although 
there is a fairly large scatter. The top left panel of Figure [3] shows the histogram of isotropic 
energy values, Eiso, calculated using Equation [7] using the spectral parameters in Tableland 
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redshift in Table [31 A large range of values for Eiso is produced because of the spread in the 
isotropic energy density at the peak Qp^, the peak energy Ep^ and the bolo metric integral 
Iboi- The top right panel of Figure [3] shows the peak energy values, Ep^, plotted against the 
characteristic energy, Ey^^- There is a tight correlation between these 2 parameters because 
of the form of the Band function. To a first approximation Ep^ = 0.23Ewz (the solid line 
in Figure [3]) although the best fit correlation is a little steeper {Ep^ oc El;l'^) and the small 
scatter evident in Figure [3] is introduced by differences in the spectral indices, Px and P^. 
In fact, the bolometric integral is well approximated by a function of the form 

hoi ^ Ifu = Elt'^ exp(co + C2(3x + c^P^) (10) 

where the coefficients Co,Ci,C2,C3 can be found by a least squares fitting procedure. A 
comparison of Ifu and hoi for the GRBs listed in Table [T] is shown in the bottom left panel 
of Figure [3] together with the best fit coefficients. We can use If a in place of hoi and estimate 
Eiso ~ Efit. The distribution of the ratio Eiso/Efu is shown in the bottom right-hand panel 
of Figure [31 For the majority of objects the estimation, Efu, is within ±10% of the value 
obtained by numerical integration. There are a few GRBs with a larger discrepancy but all 
are within ±20% which is a very small perturbation in comparison with the dynamic range 
of the Eiso values. 

Using hoi ~ I fit we can express Eiso as an explicit function of Ep^: 

E,so ^ Qp.E^^^'-^'' exp(l ± Co ± (C2 - l)Px + c^P^)- (11) 

The immediate origin of the Amati relationship is now clear. Given Equation [Til some degree 
of correlation between Epz and Eiso is guaranteed. The nature and spread of this correlation 
will depend on the relationship between the fiux density, Qpz, and the peak energy, Ep^, 
and the distribution of spectral index Px- It could be that Qpz and Ep^ are correlated in 
such a way to cancel the apparent dependence on Epz but this is highly unlikely. This 
correlation arises because the GRB spectral profile has the form of Band function (Equation 
[1]) with a particular range of values for the spectral indices, Px, P-y, and the energy Ec. So 
understanding where the Amati relation comes from is really the same as understanding why 
the spectra have this functional form in the first place. 

Figure [H shows the Amati relationship for the GRBs in Table [H Here and subsequently 
we use the exact form for Eiso, calculated from hoi, and not the approximation involving 
I fit which was only introduced to derive Equation [TTl The correlation line shown (derived 
ignoring the obvious outliers) is Ep^ cx E'^^^^ consistent with Amati 2006, Ep^ oc E'^^^. All 
the short bursts are outliers with low Eiso values compared with the long bursts of similar 
Epz value. The other notable outliers are GRB980425 and GRB060218 (see Amati 2006, 
Campana et al. 2006). The XRFs (characterised by the hardness ratio of the low energy 
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spectra, see below) all fall on the lower edge of the correlation with low Epz compared 
with Eiso- A more fundamental relationship is that between the flux density Qpz and the 
characteristic energy E^^z which is also shown in Figure IH It appears that, disregarding the 
short bursts, the Amati correlation is tighter than this new relationship but this is deceptive. 
Unlike Eiso and Epz, Qpz and E^z are independent and their product provides the isotropic 
energy Ei^o (Equation [7]). We now have a correlation which goes beyond the simple fact that 
GRB spectra have the Band function profile. Crudely, Qpz is a measure of the height of the 
spectrum as plotted in Figured] and E^z (which is itself a function of Epz, Px and P^) is a 
measure of the characteristic photon energy. There is a weak correlation between these two 
quantities, E^jz oc Qpf, as can be seen in Figures H] and [TJ However, the pattern of outliers is 
the same as for the Amati relationship. The short bursts, and the sub-luminous long burst, 
GRB980425, have significantly low Qpz values but E^z values which are comparable to the 
gamut of long bursts. The bursts designated as XRFs (see below) all lie in the low tail of 
the Ey^z range but have Qpz values which are similar to many long bursts. 



4. The rate profile and luminosity time of the prompt emission 

The analysis above has highlighted the well known problems associated with the Amati 
relation and other correlations involving Eiso. We now consider a way of converting E^go 
into a characteristic luminosity to see if this can improve the situation. The variety of 
time variablity in the prompt emission from GRBs is astonishing. Some bursts consist of 
a single Fast Rise Exponential Decay (FRED) profile, other have multiple peaks, some are 
very spikey with rapid variations while others have a smoother profile. The luminosity is 
continually varying between bright, short peaks and low troughs and in some cases the flux 
drops below the detection threshold for a while before flaring up again. With such a range 
of behaviour defining some characteristic luminosity and/or duration is tricky. Reichart 
et al. (2001) showed that the peak luminosity correlated with a variability measure V 
computed by taking the difference between the light curve and a smoothed version of the 
light curve where the smoothing or correlation time was the time taken to emit the brightest 
fraction / of the flux, T^. They showed that the most robust correlation was obtained for 
/ 0.45. The correlation of the peak luminosity with T45 has been adopted by subsequent 
authors, for example Guidorzi et al. (2005), Firmani et al. (2006), but in all cases the peak 
luminosity must be defined using some small arbitrary bin size (typically 1 second) and the 
only connection between the total fiuence and the peak luminosity is indirect, through the 
T45 value. 

The variability measure V depends on the correlation of structures (peaks, troughs etc.) 
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in the light curves. Here we try a different approach in which the sequence of features or 
events in the hght curves is abandoned completely. We identify the time periods in which 
significant flux is measured and then construct a rate profile by sorting the sequence of count 
rate samples from these time periods into descending order to produce, for every GRB, a 
monotonically decreasing function, fsits), where tg is sorted time. The total sum of all the 
samples should be the total count fiuence and the profile is normalised by dividing by this 
fluence so that the integral under the profile is unity. Such a rate profile shows what fraction 
of the burst is spent at what fraction of the peak rate and has the general form shown 
schematically in Figure [51 Examples of these rate profiles are shown in Figures [6] and [HI 
The time periods in which significant flux is detected were found by successive correlation 
with boxcar functions of increasing width. It doesn't matter if the total duration of these 
periods is a little larger than required to capture the total fiuence because the small excess 
of samples in the tail can be dropped and the rest of the profile is unchanged. Remarkably 
the shape of these rate profiles is surprisingly similar for all GRBs and is insensitive to the 
time bin size used as long as it is not too large or too small. If the bin size is too large 
then there may be too few samples defining the profile, but we found that a number of bins 
> 20 was fine. Using excessively large time bins can also hide significant real structure in 
the fluctuations of the light curve and this should be avoided. At the other extreme, if the 
bins are too small the number of counts per bin may drop to single figures and the profile 
shape is again compromised. In practice all long bursts are well represented using ^ 64 ms 
bins while short bursts require ~ 4 ms bins or something similar. 

The influence of statistical fluctuations (noise) on the rate proflles is rather strange. Be- 
cause the integral is normalised to unity statistical fluctuations on the total fluence are not 
included. The proflle reflects the distribution of the detected flux over a range of brightness 
but is not influenced by uncertainties in the total flux. The sorting of bins into decreasing 
brightness order also ensures the proflles are always smooth with the larger errors or dis- 
tortion due to noise accumulating at the start and end of the proflle. This is often most 
noticable as a slight increase in gradient or curl over at the end of the proflle. Although 
errors can be estimated for each of the samples, di, Chi-squared minimization using these 
errors cannot be employed for any function fltting because the sorting operation destroys 
the meaning of the errors, i.e. the scatter of the sorted data values about the fltted function 
is not governed directly by the errors on di. 

Most proflles are well represented by an empirical function of the form 



where Te is the total emission time or duration of the proflle, is the level of the proflle 




(12) 
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at Te and represents the minimum detectable flux (or luminosity) and Cl is a luminosity 
index which describes the curvature. This function is illustrated in Figure [51 Because the 
profile integral is normalised to unity the peak value at the start is fo + Ie = I/^l, where 
Tl is a luminosity time in seconds. The peak flux is then given by the fiuence divided by 
the luminosity time, Ftot/TL cts s~^ or, perhaps more intuitively, the peak flux multiplied 
by Tl is the total fiuence. Because Tl is derived from the functional fit of all the data it 
does not depend strongly on the time bin size (as discussed above) and therefore the peak 
flux calculated using is also not dependent on the binning. If Cl = ^ the profile is linear 
and if Cl > 1 the profile is concave and the fraction at high rate is smaller. If Cl < 1 the 
curvature would be negative but this is not seen for any GRBs. So Cl is a measure of the 
sharpness or spikiness of the profile. 

We fitted all rate profiles with the function fg given by Equation [T2] finding the best fit 
values for the parameters Tl, Cl and using a least squares statistic 



where is the total number of samples, di, of time width At. Note Te is fixed as the 
cumulative duration of all the significant samples detected, Te = NAt. The S statistic has 
properties similar to reduced Chi-Squared, with typical values in the range 0.5-2.0 (set by 
the scaling factor of 100), independent of the number of samples or the sample size At. 
Table [2] provides a complete list of all the temporal parameters. This table also includes the 
instrument and a GRB classification using the usual observational definitions: Short bursts 
if Tgo < 2 s and X-ray Flashes (XRFs) if fiuence(l-30 keV) /fiuence (30-500 keV)> 1. 

Figure [6] shows examples of typical fits. Note that sorted time is scaled by 1/Tl; and the 
fs values by Tl so that both axes take the range 0-1. The top-right panels show GRB021211 
which is a typical FRED burst and has a low curvature index, Cl = 1.28. The top- left and 
bottom right panels show GRB990510 and GRB070521 which have more complicated flaring 
structure but are well fitted with Cl values of 2.48 and 1.69 respectively. The remaining 
objects have short bright spikes and extended low level emission, a class discussed by Norris 
& Bonnell (2006). GRB050724 and GRB051221A are essentially short bursts followed by 
a low level, extended tail and the combination of these features produces large Cl values, 
3.17 and 2.85 respectively. For GRB051221A, E = 3.22 which is rather high. In this case 
the short spike followed by the extended tail produces an extra feature or wiggle in the rate 
profile which is not fitted by the simple fg function. Equation [T21 For these and similar 
bursts a sample size of 4 ms was used to accomodate the profile of the initial short spike. 
The left-hand panel of Figure [7] shows the distribution of E and Tl values for all GRBs 
in Table [21 There is no correlation between the goodness of fit measured by S and the 
luminosity time, Tl. The same is true for S and the luminosity index Cl- Figure [H shows 




(13) 
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the worst fits of rate profiles with large S values. In all these GRBs the peak value, l/T^, is 
a good approximation to the data peak but the fit is compromised by undulating features. 
GRB990705 and GRB061007 represent a small group of bursts which have fiares that rise 
fast, are reasonably fiat at the top and decay fast. These produce a characteristic S-feature 
in the profile. Only 12 rate profiles (out of 101) have E > 2 and only 3 of these have a 
substantial mis-match, GRB990705, GRB061007 and GRB061210. The latter is an extreme 
example of a short burst, Tl = 0.03 s, which has an extended low fiux tail giving Tqq = 85.3 
s. We note that GRB991216 has a faint pre-cursor just visible on the lightcurve plot. 

The combination of luminosity time, Tl, and curvature index, Cl, gives us information 
closely related to T45. The right-hand panel of Figure [7] shows the correlation between Cl 
and the ratio of T45 calculated directly from the sample values di and Tl from the fitted 
function. T45 could be calculated by integration of the fitted function using the parameters 
Tl, Cl, /e and Te and this would produce a smooth curve of Cl vs. T^^/Tl if fs were zero 
or constant. The parameter Cl, for example, could be replaced by T45 and the fitted function 
would still be uniquely defined. The scatter in Figure [7| results from the small differences 
between the data and the fitted function and the value of Je which is generally much smaller 
than 1/Tl but different for each GRB. Error ranges for Tl and Cl were estimated assuming 
the statistic E has properties similar to reduced Chi-Squared. The errors so derived are not 
statistically correct, because of the odd statistical nature of the sorted rate profile, and in 
some cases they are an over estimate as is evident from the scatter in Figure [3 

Although the minimum fiux level, Je, was included in the fitting it is a measure of the 
instrument sensitivity rather than some intrinsic property of the rate profile. If the noise 
level were lower the number of significant samples detected would increase, Te would get 
bigger and /e would decrease. The instrument would detect a slightly larger fiuence, Ftot, 
and the fitted value of Tl would increase a little, however, the peak fiux level, Ftot/TL would 
remain unchanged and Cl would be essentially the same. The analysis of the rate profile 
described above provides a robust estimate of the peak fiux (or peak luminosity) using all the 
available light curve data and is not biased by the instrument sensitivity providing the burst 
detection significance is secure in the first instance. The error on the peak fiux so estimated 
is dominated by the error on the fiuence rather than any error associated with estimating 
the luminosity time, Tl- It is also unchanged by the choice of sample size. At, providing 
the number of samples is sufficient to capture the details of the emission profile as already 
discussed above. We can never be sure that resampling a light curve with a smaller At will 
not reveal a very short, bright, isolated spike which was hidden by the previous binning and 
this would compromise the shape of the profile, but such has not been seen in any of the 
GRB light curves analysed so far (about 250 including all Swift bursts to date). 
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Using the redshift, we can calculate the luminosity time in the source frame, T^^ = 
+ ^) and 90% duration in the source frame, T902 = T^i^jiX + z). The peak luminosity 
multiplied by the Tj^^ gives the isotropic energy, LisoTiz = Eiso ergs. This simple property 
of makes it a highly significant measure of the burst duration and is why we chose to 
call it the luminosity time. Such a time is often introduced in theoretical dicussions, see for 
example tj in Thompson et al. (2007) or ti^urst in Ghirlanda et al. (2007). Above we have 
described a method to calculate this time for every GRB. 



5. Characterisation of the prompt emission in the source frame 

The prompt emission of each GRB in the source frame is characterised by the peak 
energy density, Qpz ergs keV~^, the characteristic photon energy, E^z keV (which embodies 
the spectral indices [3x, P-y and the peak energy Epz, Equation [8]), the luminosity time, 
Tlz s, and the luminosity curvature, Cl- Figure [3 shows T^z plotted against the standard 
measure of burst length Tgo^ where the dashed line shows equality. For bursts consisting of 
a single smooth pulse then Tlz ~ ^902 • If there is more structure in the light curve and, in 
particular, if there are periods when the flux drops to zero then Tlz < Tgoz- In some cases a 
short precursor pulse is followed by a long time gap before the main burst starts and then 
Tlz « TgQz- So the ratio of the two times is a crude measure of the variability but this 
includes all time scales and long periods when no flux is detected and is not equivalent to the 
short time scale variability defined by Reichart et al. (2001). The top-right panel of Figure 
[9]shows the distribution of Tlz- Two peaks containing the short-bursts, centred around 0.05 
seconds, and long-bursts centred at 5 seconds, are clearly visible. The distribution of C^, is 
shown in the lower left-hand panel of Figure M Most bursts are contained in a symmetrical 
peak centred on Cl = 1.6. The few bursts with Cl > 2.2 include the short bursts which 
have a long weak tail and bursts which exhibit several very short spikes on top of a more 
generally smooth emission. The bottom right-hand panel of Figure [9] shows the distribution 
of the Band lower energy spectral index, Px- Hard bursts have < and softer bursts 
have Px > 0. We do not show the distribution of the high energy spectral index, P^, because 
this parameter is only available for a few bursts and in most cases it was set to = 1.3 
which is the approximate average found by BATSE. 

The distributions of the remaining parameters, the characteristic photon energy, E^j^, 
and the peak energy density, Qpz, are shown at the top of the Figure [TUl E^^j^ stretches over 
two decades from 100 keV to 10000 keV. Qpz has a much larger spread with a main peak 
spanning three decades and a low energy tail covering another three. Since the product of 
the two gives us Eiso the range of isotropic energy is very large, as is evident from Figure 
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|3] and the Amati relation plotted in Figure HI The peak energy density, Qpz, is correlated 
with the luminosity time T^^ as demonstrated by the bottom left-hand panel of Figure [TOl 
Short bursts have Qpz < 10^^'^ while in general long bursts have larger Qp^ values. The two 
notable exceptions are, as before, GRB980425 and GRB060218 which are long bursts with 
very low luminosity. Five short bursts with long tails that are classified as long because their 
T90 > 2, GRB050603, GRB050724, GRB061006, GRB061210 and GRB070714B have Tl, 
of 0.22, 0.38, 0.33, 0.02 and 0.48 s respectively and these sit below the main long grouping 
along with the shorts. The XRFs tend to have lower Qp^ and lower T^^ values within the 
long burst population. The peak luminosity density of a burst is given by Qpz/Tiz ergs 
keV~^ s~^. This has a much narrower distribution than either Qp^ or T^z with a 90% range 
just over 2 decades, 3.1 x 10^^ — 7.8 x 10^^ ergs keV^^ s"^, as is clear from the histogram in 
the bottom right-hand panel of Figure [TUJ Both short and long bursts have similar values of 
peak luminosity density (the short bursts are shown as the white histogram, all bursts are 
shown in the grey histogram). 

We have used Principal Components Analysis (PCA) to investigate the scatter within 
the parameter space described above {Tlz, Cl, Qpz and E^z)- This analysis confirms that 
there is indeed a correlation between T^z and Qpz and the best fit is very close to pro- 
portionality with index 0.89 but there is considerable scatter with Pearson's correlation 
coefficient is r = 0.5, Kendall's r = 0.3, 4.5a (see Figure [TOj) . If we fix this index to 
unity then the only other significant correlation is between Qpz/T^z and the characteristic 
energy Ey^z- This is shown in the left-hand panel of Figure [TT] including the best fit correla- 
tion, Eyjz ^ (Qpz/Tiz)^''^^ which has Pearson's correlation coefficient r = 0.55 and Kendall's 
T = 0.32, significance 4.8cr. The critical difference between this plot and the right-hand panel 
of Figure H] (the Amati relation) is that the peak energy density Qpz has been converted to 
a peak luminosity density by dividing by the time Tlz- The large difference between the 
short and long bursts has disappeared and most bursts are now clustered in a small area 
on the energy-luminosity plane. It seems that all correlations involving the properties of 
GRBs must have outliers and this is no exception; GRB980425 still refuses to conform but 
the remaining 100 bursts come into line. 

The correlation shown in Figure [TT] between the peak luminosity density and the char- 
acteristic photon energy in the source frame is the first GRB relationship to unify the short 
and the long bursts. If the peak luminosity density is multiplied by Ey^z the x-axis becomes 
the peak isotropic luminosity, Liso ergs s^^. The correlation of E^jz vs. Lj^o is shown in the 
right-hand panel of Figure [TT] along with the best fit 
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which has a Pearson's correlation coefficient of r = 0.75, Kendall's r = 0.54, significance 
8.0(7. Thus Figure [11] encapsulates a major result of this work, showing a high quality 
correlation of characteristic photon energy with peak isotropic luminosity for 101 GRBs 
including 9 short bursts and 7 XRFs. The correlation between E^^ and L^^o is similar to 
those reported by Yonetoku et al. (2004) and Firmani et al. (2006) but there are important 
differences. Here we have estimated the peak isotropic luminosity from the rate profile so we 
are not restricted to long bursts or a particular time bin size, and the peak energy, Ep^, is 
replaced by the characteristic photon energy, E^^- We note that the correlation derived by 
Yonetoku et al. (2004) is significantly steeper, Ep^ cx but they used a rather small 

sample of 16 GRBs. We can identify 13 of these objects in our sample and we find they 
give Epz oc L°,^^^° °^ with Pearson's correlation coefficient r = 0.81, consistent with their 
result. The same 13 objects also give E^^ oc L^^^^^"'""^ with r = 0.88 so using E^^, in place 
of Epz gives a slightly tighter correlation with a shallower slope which is consistent with our 
result obtained from the full sample of 101 bursts. Unlike the Firmani et al. relationship 
the present correlation does not contain T45^. We tried including T45Z in the PCA but found 
no significant correlation or reduction in the scatter. If we replace Ey^^ by Ep^ in the PCA 
of the complete sample then Ep^ cx L^;^^^ with r = 0.71 and Kendall's r = 0.52, significance 
7.7(7 so, again, using E^z yields a tighter correlation with a shallower slope compared to Epz. 
The small change in slope arises because the correlation of Epz with Ey^z is not quite unity 
(see Figure [3]). 

For each burst we calculate Kz which is a measure of its displacement perpendicular 
from the the best fit correlation line in the right-hand panel of Figure [TTl 



This is a function of the ratio of the characteristic photon energy to the peak isotropic 
luminosity. The constants quoted in this definition are the mean values of the parameters 
so they represent the centre of the clustering of objects within the parameter space. The 
distribution of Kz is plotted in Figure [T2l The mean value is \ogiQ{Kz) = or, equivalently, 
Kz = I. Hard-dim bursts (including most short bursts) have Kz > I, soft-bright bursts 
(including all XRFs) have Kz < 1. The distribution is approximately log- normal (the best 
fit Gaussian profile is shown in Figure [T2!) and has a rms width of cr[log]^o(-^z)] — 0.19. 90% 
of the GRBs (90 objects) are contained in the range 0.45 < Kz < 1.95. The obvious outlier 
is GRB980425/SN1998bw which is either very sub- luminous or has an exceptionally high 
peak energy for such a dim burst. Under the hypothesis that Kz is constant, = 421 with 
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99 degrees of freedom and the mean of the estimated errors on log 10{Kz) is 0.12 so there is 
clear evidence for intrinsic scatter in with an estimated 90% range of 0.57 < < 1.75. 
The largest uncertainties arise from the estimation of E^z because this depends on Ep^ and 
the spectral indices Px, P-y- The mean value of for the pre-Swift bursts is —0.05 and 
for Swift bursts is 0.02 so they are statistically indistinquishable. The distribution for pre- 
Swift bursts, plotted as the white histogram in Figure [121 sits symmetrically within the total 
distribution. The right-hand panel of Figure [12] shows log^gl-^z) a function of redshift, 
z. There is no obvious trend. The objects with the smallest errors that contribute most to 
the high show no dependence on redshift. Table [3] provides a complete listing of the rest 
frame parameters and associated errors. 



6. Discussion 

Within the new parameterisation of the temporal and spectral properties of the prompt 
GRB emission the three important quantities are the characteristic energy in the source 
frame, Ey^z keV (Equation [8]), the energy density at the peak of the E.Fz{E) spectrum, Qpz 
ergs kev~^ (derived from the total fluence. Equations [21 [5] and [6]) and the luminosity time, 
Tlz s, derived from the rate profile. The ratio Qpz/T^z gives us the peak luminosity density 
in ergs keV~^ s~^ where "peak" corresponds to both the maximum in the E.Fz{E) spectrum 
and the maximum flux level in the light curve. Ey^z and Qpz/T^z are correlated and it is this 
correlation which gives rise to the Amati relation (Amati et al. 2002). The instantaneous 
maximum brightness of the prompt emission is characterised by a function of the photon 
energy/peak luminosity ratio Kz given in Equation [151 This is not a constant but it covers 
a remarkably small dynamic range compared with the constituent parameters, E^z, Qpz and 
Tlz- Given the measurement errors it is difficult to make an accurate estimate of the intrinsic 
dynamic range but it is certainly less than 0.5 < Kz < 2.0 and this holds for 100 GRBs in 
the sample of 101 we have analysed including long, short and XRFs, the exception being 
GRB980425. 



6.1. Is Kz intrinsic? 

We might wonder whether the narrow range in Kz is an artefact of the observational 
data or something instinsic to the nature of the GRB emission? An artificial tightness of 
the energy-luminosity correlation could arise in several ways; the observed quantities may be 
correlated by some property of the instrumentation/measurement, the measured positions 
of GRBs in the energy-luminosity plane could be incorrect because of some systematic er- 
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ror/bias or GRBs from certain areas in the plane may be selectively missed. The measured 
quantities in the observer frame which map to E^^ and Qpz/Ti^ are the peak energy, Ep 
keV, and the spectral energy density at the peak 

U = ^exp[(l - 1)]E;^- (16) 

J-L 

ergs cm~^ keV~^ s^^. These are plotted in the top left-hand panel of Figure [131 There is 
no tight clustering or significant correlation. Pearson's correlation coefficient is r = 0.27 
and Kendall's r = 0.13, significance 2.0a. The 90% range of fp is 4.1 x 10"^° - 2.8 x 10"^ 
ergs cm~^ keV~^ s~^ and the 90% range of the observed peak energy is 34 — 412 keV, both 
one to two orders of magnitude. Using the redshift to transform these into Qpz/Ti and 
E^z produces the distributions shown in Figure [101 These quantities have slightly narrower 
distributions in the source frame with 90% ranges of 3.1 x 10'^'' — 7.8 x 10^^ ergs keV^^ s^^ 
and 511 — 3450 keV respectively. Finally they combine in Kz which has a rather narrow 
90% range of 0.45 — 1.95 and some of this is attributable to the measurement errors. It is 
very unlikely that some systematic error or bias in the measured quantities which have a 
large dynamic range and are not correlated conspires to give such a tight correlation and we 
conclude that Kz encodes real, useful, information about the source frame properties of the 
prompt emission. The rate profile fitting not only provides us with the peak fiux density 
level but also the minumum detected fiux density 

fm = iVtot/i^ exp[(l - M{E;' - 1)] V (17) 

ergs cm~^ keV~^ s"^. The top right-hand panel of Figure [131 shows fm vs. fp. There is 
some clustering of the weaker bursts along the line fp ~ 5/^ and clearly the area above 
this line in the top left corner is below the threshold. It could be that we are preferentially 
missing hard-dim bursts while soft-dim bursts are detected but redshift works in our favour 
because the distant dim bursts are redshifted into the lower observation energy band where 
the sensitivity is higher and time dilation stretches the light curve so we have longer to 
detect the emission. We are undoubtedly missing low luminosity bursts especially at high 
redshift. The lower panels of Figure [T3l shows Liso and Eyjz plotted vs. redshift. The hard 
(high Eyjz) and most luminous (high Liso) sources are seen at all redshifts while the softer, 
weaker sources are only seen at low redshifts, entirely as one would expect, but as we can 
see from Figure [T2l there is no obvious difference in the distribution of Kz as a function of z. 
There is no reason to suspect that absence of dim bursts too weak to detect is biasing the 
distribution in the photon energy-luminosity plane. 

Our conclusions are somewhat different from Butler et al. (2007). The combination of 
recent Swift detections with pre-Swift results confirms the general correlation between Epz 
and Eiso (the Amati relation) but the spread is indeed large and many bursts, including all 



-16- 



the short bursts, are extreme outliers from the bulk correlation of the long bursts. Such 
a correlation is, in part, a simple consequence of the shape of a typical GRB spectrum 
(the Band function) but the spread and presence of many outliers renders this correlation 
insensitive in testing of cosmological world models. There is a clustering of events in the 
ratio of fluence/duration (effectively luminosity), however, we think the root of the problem 
is not observational bias or sensitivity thresholding but rather that the Amati relation (and 
similar correlations involving Eiso) are looking at the wrong parameter space. Yes, this is the 
demise of the existing pre-Swift high-energy correlations, but if we re-cast them in terms of 
the instantaneous peak luminosity and we replace Epz by Ey,z, which combines the spectral 
parameters Ep^, (5x and then they reappear in a new light. The quantities in the observer 
frame are not correlated, the short and dim bursts are no longer outliers in the source frame 
correlation and there is no difference between Swift and pre-Swift detections. 



6.2. Correlations involving evolution of parameters 

All the analysis presented here involves average spectral and temporal properties of 
the prompt emission. Ey,z and Qpz are derived from the time-integrated spectra and Tlz is 
estimated from the full energy band light curves. We know that GRB spectra evolve with 
time and the light curves are different in different energy bands. In general the spectra 
soften as the burst proceeds, (3x increases and Ep decreases with time (e.g. Goad et al. 
2007, Page et al. 2007). The light curves are shorter and more spikey at high energies than 
they are at low energies (Reichart et al. 2001). The lag- luminosity (Norris et al. 2000) 
and variability-luminosity (Reichart et al. 2001) correlations are testament to this temporal- 
spectral evolution. In this work we have estimated the hardness and brightness of just the 
peak emission. If instrumentation could follow the evolution of the characteristic energy and 
luminosity through the light curve each burst would form a track on the energy- luminosity 
density plane which may run from top right to bottom left with Kz ~ constant. The lag 
and variability correlations may provide a means by which scatter can be introduced in 
although physical reasons for this are not immediately apparent. Further analysis and better 
quality data are required to explore the evolution of through individual bursts. 



6.3. Emission processes 

The correlation between the hardness and brightness of GRB spectra, previously in the 
form of the Amati relation and now the correlation between 5',^,^ and Qpz/T^z, is a challenge 
to theoretical modelling of the prompt emission. Within the standard fireball picture there 
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are many variants involving internal and external shocks in which synchrotron emission, 
inverse Compton scattering and photospheric emission feature, and the fireball itself may be 
dominated by kinetic energy or magnetic energy (Poynting flux). The initial problem is to 
predict a spectrum which has the general form of the Band function with a spectral break 
or curvature characterised by some energy, Ec, Ep or E^z, and the second problem is to 
predict the coupling between the characteristic energy or hardness of the spectrum and the 
luminosity (see the review by Zhang & Meszaros 2002). 

With a kinetic energy dominated outflow and a simple synchrotron model generated by 
internal shocks, incorporating a peak in the electron energy one expects Ep^ oc r^^t~j^L^/^ 
where T is the bulk Lorentz factor, tyar is the typical variability time scale associated with 
the internal shocks and L is the luminosity (Zhang & Meszaros 2002). This is consistent 
with the Amati relation if L oc Eiso (which is not the case if we include both short and 
long bursts as shown above) and there is a constancy of both F and t^ar across all bursts, 
which seems unlikely (Rees & Meszaros 2005). The relationship E^^ Lf^^ derived here 
and shown to hold for all bursts is significantly flatter than the Amati relation. If the 
Lorentz factor depends on luminosity, V oc , then we can choose /? ~ 1/8 to match the 
observed correlation providing t^ar is independent of luminosity and approximately constant 
for all bursts. It is not obvious why the Lorentz factor should have such a specific and low 
dependence on luminosity and, again, why t^ar should be constant when the burst durations 
(Tgox or Tlz) have such a large dynamic range. Within this model the radius at which 
the emission occurs is given by r ~ ctyar^'^ so if we assume typical values of t^ar ~ 0.01 s 
and r ~ 300, r ~ 3 x 10^^ cm. Furthermore, by considering the onset of X-ray afterglows 
observed by Swift Kumar et al. (2007) estimate that emission originates at much larger 
radii, between 10^^ and 10^^ cm, and suggest that synchrotron/inverse Compton parameters 
cannot account for the prompt emission. 

Alternatively, we can can consider a thermal origin for the peak in the E.F(E) spectrum 

and the correlation of the characteristic energy with luminosity, see for example Rees & 
Meszaros (2005), Ryde (2005) and Ghirlanda et al. (2007). If the photosphere of the 
expanding fireball has radius Rq, Lorentz factor Fq, a blackbody spectral component with 
temperature T^^ and isotropic luminosity fraction Ebb of the total isotropic luminosity Lj^o, 
then the observed temperature, Tobs = (4/3)FoTfefe, is given by 

Tobs ^ ( To W cm \ / BbbL^so \ . . 
1460 keV \ Ro J VlO^^ergss-V ^ ' 

Thompson (2006). This is just the Stefan- Boltzmann law modified to account for the 
relativistic expansion rate of the photosphere and it matches the observed correlation if 
£■^2 oc Tabs and eH^Tq/Ro is approximately constant for all bursts. The observed spectrum 
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is not a single temperature blackbody. If Ebb << 1 then any single temperature blackbody 
component is diluted by non-thermal (possibly power law or inverse Compton) components 
which combine to give the Band function. If e^f, ~ 1 the observed spectrum must result from 
the summation of a large number of thermal components with a spread of temperatures that 
give an average of Tobs, (for example, in a manner similar to that described by Ruffini et al. 
2004). We can re-arrange Equation [18] as 



If Tobs ~ Ewz this is essentially the same as the definition of the photon energy- luminosity 
quasi- const ant, Kz, and so, under this interpretation, the scatter in the correlation arises from 
variations in the fireball dimension, Rq, the Lorentz factor, Fq, or the blackbody luminosity 
fraction. Ebb- If = 1 then we have an average fireball with Rq/{eI[^Tq) ^ 10^ cm. 

1 /2 

If Kz > 1 the fireball has a higher than average E^^l Tq and/or a smaller radius, Rq. If 

1 /2 

Kz < 1 the fireball has a low ^j,^ Fq and/or large radius. With Fq = 300 and Ebb = 1 then 
i?o ~ 3 X 10^ cm which is the thermalization radius (i.e. the radius of the jet or fireball 
photosphere) estimated by Thompson, Meszaros & Rees (2007). This radius is much smaller 
than estimates arising from the internal shock model or recent estimates involving the onset 
of the X-ray afterglow (see above) so independent estimates of the radius of the prompt 
emission and/or the Lorentz factor may help to discriminate between thermal and internal 
shock models. We also note that baryonic photospheres are governed by physical argument 
(Meszaros et al. 2002) such that the ratio Rq/^o is constrained and the relationship between 

1 /2 

Ef^l Tq/Rq and Kz may not be as simple as we have indicated. 



The Eujz — Qpz/Tiz correlation should be useful as a pseudo redshift indicator but there 
is evidence for intrinsic scatter in Kz, the errors in determining the characteristic energy and 
luminosity density are large and the correlation gradient is rather low, 0.24. In particular, 
reliable estimation of Eu,z requires an accurate measurement of the broadband spectrum 
and always involves some extrapolation to cover the source frame energy band 1-10000 keV. 
Without a good measurement of the peak energy, Ep, we can't calculate a good estimate of 
Qpz from Qz- The launch of GL A ST in the near future will hopefully provide excellent high- 
energy spectral measurements which will tie down the spectral parameters more precisely. 
The intrinsic scatter may arise in several ways. If the dominant emission mechanism is 
non-thermal then the coupling of the Lorentz factor of the expansion with the luminosity 
and the variability time associated with the internal shocks may be the root cause. If 
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6.4. Kz as a cosmological probe 
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thermal processes dominate then the ratio Rq/Tq may vary as discussed above. Because 
we can identify classes which fall predominately at > 1 (shorts) and < 1 (XRFs) 
there is some hope that additional parameters which distinguish these classes may serve 
to narrow the distribution. In the above discussion we made no mention of coUimation or 
beaming of the outflow. Since the Amati relation involving Eiso has been transformed into the 
present universal correlation involving Liso the simple beaming argument that underpinned 
the Ghirlanda relation (Ghirlanda et al. 2004) is not directly applicable and currently there 
is no simple physical model which links collimation to scatter in peak luminosity or the 
characteristic photon energy. However, it is not unreasonable to suppose that collimation 
may introduce scatter in the peak luminosity, and that correlation of with afterglow 
parameters such as optical jet break times tbreak or the time of the start of the final X-ray 
afterglow, (Willingale et al. 2007), may be fruitful. In its present form is not a 
sensitive cosmological probe but the signs are that it may be in the future. 



7. Conclusion 

The equivalent isotropic energy, Eiso ergs, of a GRB can be expressed as the product of 
two source frame terms, a characteristic photon energy, E^^ keV, calculated from the shape 
of the spectrum across the range 1-10000 keV and the energy density at the peak of the 
E.Fz{E) spectrum, Q^^ ergs keV~^. The correlation trend between Ey^^ and Qpz gives rise to 
the Amati relation. By stacking the samples of a GRB light curve into descending order we 
can construct a rate profile. The functional form of such rate profiles is common to the vast 
majority of bursts. Fitting the profile gives us a luminosity time, T^^ s, a measure of the 
burst duration which can be used to convert the energy density at the peak to a luminosity 
density at peak, Qpzl^Lz ergs keV^^ s^^ We can calculate the peak equivalent isotropic 
luminosity as a product Li^o = E^^Qpz/Tiz = Eiso/Tiz ergs s"^ 

Eyjz is a characteristic photon energy or a measure of the colour or hardness of the 
burst and Qpz/T^z is a measure of the instantaneous peak brightness. We have gathered 
and analysed sufficient spectral and temporal data from 101 bursts to produce the relation 
between Eyjz vs. QpzjTiz and E^z vs. Lj^o, shown in Figure [TT| which constitutes the closest 
thing we have to an intrinsic colour-magnitude diagram for the peak emission from GRBs, 
Eyjz L^sl^. All bursts are clustered such that we can construct a intrinsic colour-magnitude 
quasi constant Kz, which is a function of the source frame characteristic photon energy/peak 
luminosity ratio given by Equation [151 The range of equivalent isotropic energy that drives 
the expanding fireball is very large, 6 orders of magnitude (Figure [3]), but the instantaneous 
hardness/brightness of the peak emission covers a very small intrinsic dynamic range, ~ 4. 
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The existence and form of indicates that the physical mechanism for the Gamma-ray 
production at the photosphere of the fireball is common to all bursts and is probably thermal 
although many other possibilities are not ruled out. If the prompt spectra are dominated 
by thermal photons the scatter in may be attributed to variations in the size and/or 
Lorentz factor of the fireball. XRFs have low Tq and/or large radii. Short bursts have high 
To and/or small radii. The relation between Tl^ vs. Qpz clearly separates short from long, 
but both classes have the same instantaneous peak hardness/brightness. 

We gratefully acknowledge funding for Swift at the University of Leicester by STFC. We 
thank the authors of the BATSE {cossc.gsfc.nasa.gov/docs/cgro/batse), HETE {space.mit.edu/HETE) 
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photon energy source frame E keV 



Fig. 1. — The 1 keV to 10 MeV source frame spectra of GRBs listed in Table [U The observed 
energy band is shown as the solid line in each case. Qpz values are marked; solid dots for 
long GRBs, solid stars for short GRBs and solid triangles for XRFs. 
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photon energy source frame E keV 

Fig. 2.— The 1 keV to 10 MeV E.F,{E) source frame spectra of GRBs listed in Table [H 
The observed band is shown as the solid line in each case. Ep^Qpz values are marked at the 
peak energy Epz] solid dots for long GRBs, solid stars for short GRBs and solid triangles for 
XRFs. 
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Fig. 3. — Top panels: The distribution of Eiso and the correlation of Epz vs. E^z for GRBs 
listed in Table [U solid dots for long GRBs, solid stars for short GRBs, solid triangles for 
XRFs. The solid line is Ep^ = 0.23E^z rather than the best fit correlation which is a little 
steeper (see text). Bottom panels: Comparison of the functional fit to the bolometric integral 
I fit and the value calculated by numerical integration of the Band function over the interval 
1-10000 keV for the GRBs in Tabled! The right-hand panel shows the distribution of the 
ratio of Eiso to the value obtained using I fit. 
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'°9io(Eiso srgs) l°9io(Qpz ergs keV ^) 



Fig. 4. — Left-hand panel: The Amati relation for the GRBs in Table [TJ Right-hand panel: 
Ewz vs. Qpz for the same GRBs. Solid dots for long GRBs, solid stars for short GRBs, solid 
triangles for XRFs. 




Fig. 5. — The rate profile function. Equation [T2j 
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Fig. 6. — Typical prompt emission light curves and the corresponding rate profiles. The 
light curves have been smoothed with a boxcar function of width T45 for display purposes. 
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Fig. 7. — Left-hand panel: The distribution of luminosity times, Tl s, and the S statistic 
from the fit. Right-hand panel: The correlation between rate profile index Cl and the ratio 
of T45 derived directly from the data and derived from the profile function fit. 
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Fig. 8. — Light curves and the corresponding rate profiles for which the fit statistic E is 
high. The hght curves have been smoothed with a boxcar function of width T45 for display 
purposes. 
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Fig. 9. — Top panels: The correlation between luminosity time T^z and Tgo^ and the dis- 
tribution of Tlz- The dashed line indicates the equality Tlz = Tqqz- Bottom panels: The 
distributions of Cl and Px for the GRBs in Tabled! 
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Fig. 10. — Top panels: The distributions of E^z keV and Qpz ergs keV~^. Bottom panels: 
The correlation of T^z s vs. Qpz and the distribution of the peak luminosity density, Qpz/T^z 
ergs keV~^ s~^. The distribution for short bursts is shown as the white histogram. 
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Fig. 11. — Left-hand panel: The correlation between characteristic energy E^z keV and the 
peak luminosity density Qpz/Tlz ergs keV~^ s~^. Right-hand panel: E^z keV vs. the peak 
isotropic luminosity, Lj^o = E^^Qpz/Tlz ergs s~^. In both plots solid dots for long GRBs, 
solid stars for short GRBs, solid triangles for XRFs. The object with the largest E^^,, 8600 
keV, is GRB050904. 





Fig. 12. — The distribution of the photon energy/peak luminosity ratio corresponding 
to the scatter about the correlation line in the right-hand panel of Figure [Til The white 
histogram shows the distribution of the 26 pre-Swift bursts. The curve is the best fit Gaussian 
distribution. The right-hand panel shows vs. redshift z; solid dots for long GRBs, solid 
stars for short GRBs, solid triangles for XRFs. The horizontal dashed lines indicate the 
90% range. The objects with unusually high Kz are GRB980425/SN1998bw at low z and 
GRB050904 at high z. 
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Fig. 13. — Top panels: Observer frame parameters Ep and minimum observed flux density, 
fm, plotted against the peak flux density, fp. The solid line represents fp = fm- Bottom 
panels: Source frame parameters, peak luminosity Lj^o and characteristic energy E^j^ plotted 
vs. redshift z; solid dots for long GRBs, solid stars for short GRBs, solid triangles for XRFs. 
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Table 1. Observed spectral parameters; (3x^ low spectral index, Ep, peak energy keV (if 

no error quoted then calculated assuming fixed Ecut^ -^p = (1 ~ l^x) * 150 keV), (3^, high 
spectral index (if no error quoted then fixed at average of 1.3), Ftot-, fluence ergs cm^^, Ei 

to -^2? observed energy band keV. References: (0) Swift BAT ibid., (1) Golenetskii et al., 
2005, (2) Villasenor et al., 2005, (3) Barthelmy et al. 2005, Amati et al. 2006, Campana et 
al. 2006, (4) Golenetskii et al., 2005, (5) Golenetskii et al. 2005, Burrows D.N. et al. 2006, 

(6) Romano et al., 2006, (7) Amati et al. 2007, Montanari & Pian 2007, (8) Amati et al. 
2007, Golenetskii et al. 2006, (9) Page et al. 2007, (10) Ghirlanda et al. 2004, (11) 

Friedman k Bloom 2005, (12) Sakamoto et al. 2005, (13) Atteia et al. 2005, (14) Piro et 
al. 2005, Amati 2007, (15) Galassi et al. 2004, Amati 2007, (16) Schaefer 2007, Firmani et 

al. 2006. (17) Golenetskii et al. 2007. 
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Table 2. Observed temporal parameters; T90 sees, luminosity time sees, luminosity 
index Cl , luminosity profile fitting statistic S. Errors were estimated by assuming E was 

distributed as x^. 
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Long 


76.8 


27.3 ±0.9 


1 


.60 


± 


0, 


.06 


0, 


.32 


030323 


HETE-2 


Long 


19.6 


2.9 ±0.8 


1 


.36 


± 


0, 


.28 


0, 


.20 


030328 


HETE-2 


Long 


140 


45.8 ±0.7 


1 


.66 


± 


0, 


.02 


0, 


.27 


030429 


HETE-2 


Long 


24.6 


5.8 ±0.7 


2, 


.13 


± 


0, 


,19 


0, 


,94 


040924 


HETE-2 


Long 


5.0 


1.13 ±0.07 


2, 


.05 


± 


0, 


,08 


3, 


,40 


041006 


HETE-2 


Long 


24.6 


13.1 ±0.4 


1, 


.37 


± 


0, 


,04 


1, 


,63 


050126 


Swift 


Long 


25.7 


7.7 ±0.4 


1, 


.65 


± 


0, 


,09 


0, 


,39 


050315 


Swift 


XRF 


96.0 


21.6 ±0.5 


1, 


.62 


± 


0, 


,04 


0, 


,14 


050318 


Swift 


Long 


31.3 


5.3 ±0.3 


1 


.40 


± 


0, 


.09 


1, 


.11 


050319 


Swift 


Long 


149 


4.0 ±0.3 


1 


.93 


± 


0, 


.11 


0, 


.68 


050401 


Swift 


Long 


33.3 


6.1 ±0.2 


1 


.78 


± 


0, 


.05 


0, 


.19 


050505 


Swift 


Long 


63.0 


11.4 ±0.4 


1 


.50 


± 


0, 


.08 


0, 


,19 


050509B 


Swift 


Short 


0.07 


0.043 ± 0.005 


1 


.61 


± 


0, 


.20 


0, 


.47 
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Table 2 — Continued 



GRB 


instr 


class 




Tl 










s 




Swift 

O WIX t 


T ,C\W (T 
XjtJli^ 


8 8 

O.O 


Q + 1 

O.c) -J— U. X 


1 

X , 


.58 


+ 0, 


0*^ 


1 

X , 


45 

.t:U 


uuuuuo 


Swift 

O WIX t 


T ,C\T\ cr 

XjUII^ 


1 3 


f) Q4 _|_ f) no 




^ . 


.74 


+ 0, 


03 


n 


38 




HF,TF,-2 

LLUJ J- UJ £j 


Short 

kJllUX u 


22 







.26 


+ 0, 


08 

. U(J 


q 


19 

. X U 




Swift 

O WIX t 


yvxLX 


1 59 


f) 40 _|_ f) no 


q 


.17 


+ 0, 


06 

.uu 


n 


39 




kJ VV 1 1 Lj 


T lOTl O" 
XJW11&, 


1 55 

X ^J^J 


95 8 + 8 


1 

X 


.81 


+ 0, 


06 

. wu 


n 


95 




Swift 

kJ W 11 1 


T lOTi cr 

XJVjlig, 


SO Q 


6 1 + 


1 

X . 


.59 


+ 0, 


OQ 

. UC/ 


n 

U . 


99 


\JO\JO\JO 


Swift 

kJ W 11 Lj 


T ,c\w cr 

XJClig, 


8Q 


14 Q _(_ n c; 


1 

X 


.85 


+ 0, 


05 


n 

u . 


1 9 


\JO\JO ± o 


Swift 

kJ W 11 Lj 


^ n cw\^ 

kJliUi Lj 


45 


n nvQ _)_ n nnq 


1 

X 


.71 


+ 0, 




n 

u . 


4Q 




Swift 

O WIX t 


T ,c\w cr 

XJtJli^ 


Xt:t: 


11 1 -)_ n « 

X X . X _i_ yj . u 


1 

X , 


.74 


+ 0, 


1 1 

. X X 


n 

u . 


1 9 

. X Zj 




Swift 

O WIX t 


T ,C\T\ cr 

XJUllgi 


940 


1 1 « _|_ n 


1 

X . 


.75 


+ 0, 


03 

.uo 


n 

yj ■ 


1 6 

. XU 




Swift 

O WIX t 


T .cwv cr 

XJUllgi 


X / O 


57 + 1 


X , 


.62 


+ 0, 


04 


n 

yj ■ 


98 


UUUc/UO 


Swift 

O WIX t 


T ,c\T\ cr 

XjVJII^ 


90 3 


c c _|_ n q 


1 

X , 


.45 


+ 0, 


09 


n 

yj ■ 


1 5 

. XU 




Swift 

O WIX t 


T ,c\T\ cr 

XjUII^ 


4 1 


1 81 + n 07 

X .ox ziz u.u / 


1 

X , 


.65 


+ 0, 


05 

.UU 


n 

yj ■ 


89 


(Til 099 


XYtJllLio 


T ,r\w cr 

XJtJli^ 


1 Q7 

X U I 


in q _|_ n Q 

XU . C/ _1_ KJ.O 


1 

X , 


.81 


+ 0, 


04 


4 


Q1 

. C/ X 


n^l 1 DQA 


Swift 

kJ W 11 Lj 


T lOTi cr 

XJClig, 


■^7 


q _(_ n 


9 


.00 


+ 0, 


. X(J 


n 

u . 


98 




Swift 

kJ W 11 1 


T ,c\x\ cr 
XjCII^ 


1 5 


1 q _l_ n 9 

X . cf _i_ yj . 


1 

X , 


.79 


+ 0, 


1 6 

. XU 


n 

u . 


.0 X 


WtJ X X X X 


Swift 

kJ Wll L 


T ,cw\ cr 
XjUII^ 


49 6 


1 Q 4 _|_ n Q 


1 

X , 


.66 


+ 0, 


04 


n 

u . 


1 1 

. X X 


n^l 991 A 


Swift 

O WIX t 


OliUl t 


1 40 


n 1 1 1 Q nnq 

U.XXO _I_ u.uuo 


9 

Z: . 


.85 


+ 0, 


0*^ 

.UO 


q 


99 


yjyjyj ±\ju 


Swift 

kJ W IX u 


T lOTi cr 


14 4 

Xrt .rt 


4 4 + 03 

'-±.'-± _l_ U.(J 


1 

X < 


.46 


+ 0, 


1 1 

. X X 


n 

u ■ 


22 


Dfini 1 

uuux xu 


Swift 

O WIX t 


T ,c\T\ cr 

XjUII^ 


141 

XrtX 


17 c _|_ n « 

X / .U ZLl U.U 


1 

X , 


.78 


+ 0, 


06 

.UU 


n 

yj ■ 


35 

.OU 


fifini 1 fi 


Swift 

O WIX t 


T ,c\T\ cr 

XjUII^ 


1 05 

xuu 


C/. X ZIZ U.U 


1 

X < 


.57 


+ 0, 


1 9 

. X 


n 

yj ■ 


1 5 

. XU 


nfini94 

W\JW X Zjt: 


kJ VV 1 1 Lj 


T lOn cr 


321 


4 1+03 

rt. X _i_ w . 


1 

X , 


.58 


+ 0, 


1 3 

. X kJ 


n 


41 

X 


060206 


Swift 


Long 


7.7 


3.3 + 0.1 


1, 


.42 


+ 0, 


.05 


0, 


.22 


060210 


Swift 


Long 


220 


24.9 + 0.6 


1, 


.92 


+ 0, 


.04 


0, 


.61 


060218 


Swift 


Long 


2100 


20.4 + 0.9 


1, 


.82 


+ 0, 


.09 


0, 


.12 


060223A 


Swift 


Long 


11.4 


4.7 + 0.3 


1, 


.36 


+ 0, 


.10 


0, 


.14 


060418 


Swift 


Long 


95.8 


16.5 + 0.3 


1 


.93 


+ 0, 


.02 


1 


.05 


060502A 


Swift 


Long 


32.0 


12.0 + 0.4 


1 


.38 


+ 0, 


.06 


0, 


.77 


060502B 


Swift 


Short 


0.13 


0.058 + 0.006 


2, 


.41 


+ 0, 


.13 


0, 


.60 


060510B 


Swift 


Long 


276 


57. + 1. 


1 


.65 


+ 0, 


.04 


0, 


.24 


060522 


Swift 


Long 


69.3 


9.1 + 0.4 


1 


.77 


+ 0, 


.09 


0, 


.10 
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Table 2 — Continued 



GRB 


instr 


class 


Tgo 


Tl 










E 




Swift 


yvxLX 




« 7 _|_ n Q 

U. 1 _I_ yJ.O 


9 


.39 


± 


n 

u . 


07 

.U 1 


79 




Swift 


T ,C\T\ Cf 
XJiJilp, 


7 8 


9 q _|_ A o 


1 

X , 


.49 


± 


n 

yj ■ 


94 


14 

U . Xrt 


uuuuuu 


Swift 


T ,rwi Cf 

XjVJllej 


14 8 

Xrt .O 


7 q _l_ f) 4 


1 

X . 


.46 


± 


n 

yj ■ 


09 


1 5 

U. XU 


nfiflfiDTA 


Swift 

kJ W li L 


T iHTi cr 

XJUllc, 


100 

X yjyj 


1 5 2 + 4 


1 

X . 


.88 


± 


n 

u ■ 


04 


1 8 

U . X (J 


DfiDfiM 

WVJWVJ Xt: 


kJ VV 1 1 Lj 


T lOn 0" 

XJWllg, 


1 08 

X WO 


9fi 7 + 4 


2 


.13 


± 


n 

u . 


09 


4 70 


060707 


Swift 

kJ W 11 Lj 


T lOTi cr 

XjCII^ 


fifi 9 


1 n n _(_ n c; 

XU. U —1— u. o 


1 

X , 


.72 


± 


n 

u . 


OQ 


1 1 

U . X X 


0fi071 4 


Swift 

kJ W 11 Lj 


T lOTi cr 
XjCII^ 


1 1 ^ 

± ±o 


OQ n _(_ n « 

U —1— u. u 


1 

X , 


.62 


± 


n 

u . 


0^ 
.uu 


90 

U .Z(U 


0fi079Q 


Swift 

kJ W 11 Lj 


T cr 
XjCII^ 


1 1 9 


1 7 o _(_ n p; 

X I . _1_ u. o 


1 

X , 


.59 


± 


n 

U . 


0^ 
.uu 


44 

U .tit: 


OfiOSOl 

UUUOU J- 


Swift 

O Wli t 


OlitJl t 


4Q 


n on _)_ n ni 

U.ZjU u.ux 


1 

X , 


.70 


± 


n 

u . 


. X X 


51 

U.UX 




Swift 

O Wli t 


T .rwi Cf 

XjUII^ 


1 98 


oc c _|_ f) q 


9 


.17 


± 


n 

yj ■ 


01 

.ux 


1 95 


nfinQn4R 


Swift 

O Wli t 


T ,C\T\ Cf 


1 QO 


4 q _l_ f) o 


9 


.17 


± 


n 

yj ■ 


04 


53 
u.uo 


OfiOPOfi 

UUUc/UU 


Swift 

O Wli t 


yvxLX 


43 5 


lo 1 -1-0 4 

X.^.X 1 U.rt 


1 

X , 


.67 


± 


n 

yj ■ 


07 

.u / 


61 

U.UX 




Swift 

O Wli t 


T ,r\r\ Cf 


1Q 

Xc/ .U 


Q 1 + 3 

t?. X ZIZ u.o 


1 

X , 


.09 


± 


n 

yj ■ 


Ofi 

.uu 


31 

u.ox 


OfiOQI 9 


Swift 

O W 11 Lj 


T cr 


fi 
u .u 


1 CO -L n n« 

J..OO ZL- U.UU 


9 


.02 


± 


n 

u . 


05 

.UU 


1 57 

X .U 1 


0fi0Q97 


Swift 

kJ W 11 L 


1 lOTi cr 
XJUlig, 


1 
xu .u 


Q + 1 

O.iJ U. X 


1 

X , 


.92 


± 


n 

u . 


05 

.UU 


44 

U . TiT: 


Ofil 004 


Swift 

kJ W 11 L 


T lOTi cr 
XjCII^ 


fi 9 


9 4 + 01 

Zj.t: _I_ U. X 


1 

X , 


.89 


± 


n 

u . 


07 

.U 1 


98 


Ofil OOfi 


Swift 

kJ W 11 Lj 


T ,f^Y~\ cr 
XjUII^ 


1 9Q 


n 47 _(_ n ni 

U.t: 1 —1— U.UX 


q 


.12 


± 


n 

U . 


O'^ 
.uu 


1 9fi 

X . Z(U 


Ofil 007 


Swift 

O Wli t 


T 1 O' 
XJtJlig, 


74 9 


OQ o _(_ n o 

Zj^.Zj _J_ U.^ 


1 

X , 


.90 


± 


n 

u . 


01 

.U X 


5 04 

U .Ut: 


Ofil 1 1 0A 


Swift 

kJ W 11 u 


T ,nTi cr 

XJUllc, 


44 fi 

'-±'-± . W 


1 Q 4 -L n 5 

XO.rt _l_ U.(J 


1 

X . 


.57 


± 




u ■ 


07 

.U 1 


33 

U.UU 


Ofil 191 

L/WX X^ X 


Swift 

O Wli t 


T ,r\vi Cf 

XJlJllp, 


81 3 

OX .o 


c 04 _|_ ri AC 

U.Ort U.UU 


9 


.58 


± 


n 

yj ■ 


01 

.ux 


1 98 


Ofil 201 


Swift 

kJ W 11 u 


Short 

kJllUl u 


7fi 

yj. 1 yj 


n 00 + 01 

yj.^^ 1 u.ux 


1 

X < 


.96 


± 




u ■ 


07 

.U i 


55 

U.UU 


Ofil 91 

W\J X ^ X w 


kJ VV 1 1 Lj 


T jOTI cr 

XJVjlig, 


8'i 3 


030 + 001 

U.UOU _l_ u.uux 


4 


,29 


± 


n 

u . 


04 

. Ut: 


1 3 4 

X U . t: 


061217 


Swift 


Short 


0.21 


0.15 ±0.01 


1, 


.76 


± 


0, 


.16 


0.08 


061222B 


Swift 


XRF 


40.0 


10.7 ±0.5 


1, 


.70 


± 


0, 


.08 


0.43 


070110 


Swift 


Long 


87.3 


18.1 ±0.6 


1, 


.63 


± 


0, 


.07 


0.25 


070208 


Swift 


Long 


47.7 


2.9 ±0.3 


1, 


.69 


± 


0, 


.17 


1.12 


070318 


Swift 


Long 


119 


12.7 ±0.3 


1 


.91 


± 


0, 


.04 


0.14 


070411 


Swift 


Long 


116 


27.0 ± 0.8 


1 


.66 


± 


0, 


.05 


0.26 


070506 


Swift 


Long 


4.3 


2.1 ±0.1 


1 


.39 


± 


0, 


.12 


0.55 


070508 


Swift 


Long 


20.8 


6.59 ±0.09 


2, 


.07 


± 


0, 


.02 


1.30 


070521 


Swift 


Long 


38.4 


12.2 ±0.2 


1 


.69 


± 


0, 


.03 


0.42 
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Table 2 — Continued 



GRB 


instr 


class 








y; 


070529 


Swift 


Long 


108 


8.5 ±0.4 


1.90 ±0.09 


0.29 


070611 


Swift 


Long 


13.2 


2.6 ±0.2 


1.67 ±0.18 


0.28 


070612A 


Swift 


Long 


359 


32.4 ±0.9 


1.91 ±0.05 


0.27 


070714B 


Swift 


Long 


63.9 


0.48 ± 0.02 


2.87 ±0.07 


0.55 


070721B 


Swift 


Long 


345 


12.5 ±0.4 


1.87 ±0.05 


0.40 


070724A 


Swift 


Short 


0.40 


0.11 ±0.01 


1.75 ±0.16 


0.66 


070802 


Swift 


Long 


16.9 


3.5 ±0.3 


1.74 ±0.15 


0.27 


070810A 


Swift 


XRF 


10.9 


4.1 ±0.2 


1.18 ±0.10 


0.46 
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Table 3. Source frame parameters; z redshift, values from tabulations in Amati (2006), 
Ghirlanda et al. (2004) and http://swift.gsfc.nasa.gov/docs/swift/archive/grb_table/. 
References for all the Swift redshifts are provided on this WWW data table. Eyjz keV, 
ergs keV~^, T^^ sees, photon energy/peak luminosity ratio K^. 



GRB 




z 


E 






Qpz 




LiZ 




z 




970228 





695 


927 ± 320 


(0 


3 


±0 


i; 


10™ 


5.7 ±0.3 


1 


21 ±0 


96 


970508 





835 


906 ± 279 


(0 


9 


±0 


3; 


1049 


2.30 ±0.06 


1 


16 ±0 


83 


971214 


3 


420 


1981 ±431 


(1 


5 


±0 


4; 


10™ 


2.41 ±0.04 


1 


00 ±0 


51 


980425 





009 


550 ± 124 


(1 


7 


±0 


5; 


^045 


16. ± 1. 





18 ±0 


10 


980613 


1 


096 


807 ±382 


(0 


7 


±0 


4; 


1049 


5.7 ±0.7 


1 


02 ± 1 


12 


980703 





966 


2061 ±460 


(3 


4 


±0 


9; 


1049 


10.2 ±0.2 





50 ±0 


26 


990123 


1 


600 


5736 ± 732 


(7 





± 1 


4; 


10™ 


8.72 ±0.09 





37 ±0 


11 


990506 


1 


307 


3000 ± 674 


(3 


7 


± 1 


1; 


10™ 


6.07 ±0.04 





66 ±0 


35 


990510 


1 


619 


1539 ±215 


(1 


3 


±0 


2] 


10™ 


2.00 ±0.02 


1 


29 ±0 


42 


990705 





843 


1484 ± 189 


(1 


7 


±0 


4; 


10™ 


14.5 ±0.5 





91 ±0 


29 


990712 





430 


605 ± 119 


(1 


3 


±0 


7; 


1049 


2.9 ±0.2 


1 


78 ±0 


91 


991216 


1 


020 


2732 ± 614 


(2 


4 


±0 


7] 


10™ 


1.92 ±0.02 





85 ±0 


45 


010921 





450 


781 ± 173 


(1 


2 


±0 


4; 


1049 


8.5 ±0.3 


1 


08 ±0 


56 


011121 





360 


1395 ±217 


(4 


3 


±0 


9; 


1049 


13.9 ±0.5 





72 ±0 


27 


011211 


2 


140 


733 ± 123 


(1 


4 


±0 


3; 


10™ 


43. ±2. 


1 


38 ±0 


54 


021004 


2 


335 


1072 ±481 


(0 


4 


±0 


2; 


10™ 


3.4 ±0.2 


1 


27 ± 1 


32 


021211 


1 


010 


367 ± 79 


(3 


1 


±0 


8; 


1049 


1.0 ±0.1 


4 


43 ±2 


24 


030115A 


2 


500 


1368 ± 625 


(0 


3 


±0 


2] 


10™ 


2.4 ±0.1 


1 


04 ±1 


10 


030226 


1 


980 


1249 ± 283 


(0 


9 


±0 


2] 


10™ 


9.2 ±0.3 


1 


04 ±0 


55 


030323 


3 


370 


891 ± 445 


(0 


4 


±0 


3; 


10™ 


0.7 ±0.2 


2 


17±2 


50 


030328 


1 


520 


1106 ±247 


(2 


6 


±0 


7; 


10™ 


18.2 ±0.3 


1 


27 ±0 


66 


030429 


2 


650 


555 ±168 


(0 


3 


±0 


r 


10™ 


1.6 ±0.2 


2 


67 ± 1 


90 


040924 





859 


995 ± 137 


(0 


9 


±0 


r 


1049 


0.61 ±0.04 


1 


42 ±0 


46 


041006 





716 


524 ± 120 


(0 


4 


±0 


1; 


10™ 


7.6 ±0.2 


2 


11 ± 1 


12 


050126 


1 


290 


1444 ± 237 


(0 


7 


±0 


2] 


1049 


3.3 ±0.2 





64 ±0 


25 


050315 


1 


949 


652 ± 175 


(1 


3 


±0 


4; 


10™ 


7.3 ±0.2 


2 


25 ±1 


40 


050318 


1 


440 


547 ±118 


(2 


5 


±0 


7; 


1049 


2.2 ±0.1 


2 


40 ±1 


22 


050319 


3 


240 


1221 ±552 


(0 


6 


±0 


3; 


10™ 


0.94 ±0.06 


1 


63 ±1 


71 


050401 


2 


900 


2008 ±315 


(1 


9 


±0 


3; 


10™ 


1.56 ±0.05 


1 


14 ±0 


41 


050505 


4 


270 


2002 ± 200 


(1 





±0 


2; 


10™ 


2.16 ±0.09 





92 ±0 


23 


050509B 





225 


738 ± 220 


(0 


4 


±0 


2; 


1046 


0.035 ± 0.004 





65 ±0 


46 
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Table 3 — Continued 



GRB 


z 








Tlz 






n 606 




1 u. z _i_ u. 


in49 
xu 


4« _|_ 07 

Zj.HiU —1— U .U 1 


9 qi 7c; 

Zj. C/ X _1_ U. 1 u 


uuuuuo 


9 891 


41 oq _|_ coo 


{■\ _l_ f) 0^ 


xu 


91 q + n 005 

U.ZiXc/ u.uuu 


7q + n 93 

U. / iJ ZL- U.ZO 


fi'ifiTriQ 


1 60 
u. xuu 




(1 f) _1_ f) 0^ 

I X . L/ 1 yj 


10^^ 

xu 


048 + 003 

U.UrtO yJ.yJyJO 


9 qq _|_ 86 
z.oo u.ou 


n "10794 


n 957 


ooq -1- 1 14 

£j£jiJ _l_ X Xrt 


ff) _|_ f) 0^ 


-|a49 

XU 


qs -1- n 09 


rt. / Z 1 U.rtU 




Q7n 






XU 


K, 9 _|_ 9 
. Z _l_ u . z 


c;s + 91 

U.UO —1— u.zx 




1 71 n 

X . 1 xu 


140Q + oiq 


fo « _|_ n c;^ 


XU 


9 S ± 1 

z . _l_ U . X 


qc; + ^^4 


yJOyJOyJO 


f) 499 


040 + 1 1 n 

Ot:Z _I_ X XW 


I 0. X _I_ \J .0 


in48 

XU 


1 n 1 -1- n q 

X U. X _l_ U.O 


71 -(- 99 

U. 1 X _l_ u.zz 




1 800 


1 «Q1 _(_ 4«1 


« _)_ n 0^ 

I U. U _1_ u. z 
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